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]Cl ■ ABSTRACT 

Population III stars are now believed to contribute to the observed Near 

(N ■ 

' Infrared Background (NIRB) and heavy element pollution of the intergalactic 

, medium. Here we show that a Pop III contribution to the primordial ^He 

■ abundance consistent with NIRB constraints, might mask a lower value for 



the Big Bang Nucleosynthesis (BBN) He abundance. 
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^ ' As observations start to explore cosmic epochs close to the formation of the first, so-called 
Population III (Pop III) stars, theoretical models provide increasingly detailed predictions 
on the properties and evolution of these objects. Bond, Carr & Arnett (1983) and Carr 
(1994), discussed a possible contribution of Pop III stars to the primordial ^He mass frac- 
tion (Yp). The observational data at that time did not allow to check the validity of their 
hypothesis. Current experiments are able to determine Yp with a precision < 1% allowing 
further investigations of the problem. Such experiments are based on high signal-to-noise 
ratio spectra of a class of star-forming galaxies called Blue Compact Dwarf (BCD) galax- 
ies. These are low-luminosity {Mb > —18) systems undergoing an intense burst of star 
formation occurring in a very compact region (less than 1 kpc) characterized by a blue 
color and a HII region-like emission-line optical spectrum, which dominates the galactic 
light. BCDs are among the most metal-deficient gas-rich galaxies known. Their gas has not 
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been processed through many generation of stars, and thus best approximates the pristine 
primordial composition. Thus, Yp can be derived accurately with only a small correction 
for helium produced in stars. Moreover, the theory of nebular emission is understood well 
enough to convert ^He emission-line strengths into abundances with the desired accuracy. 
Yp is generally determined (Peimbert & Torres-Peimbert 1974) by linear extrapolation of 
the Y— 0/H correlation to O/H=0, where Y and 0/H, are respectively the ^He mass frac- 
tion and the oxygen abundance relative to hydrogen of a sample of dwarf irregular and 
BCD galaxies. Based on a relatively large sample of 45 BCDs, Izotov & Thuan (1998) found 
YpjT = 0.2443 ± 0.0015 and a scaling dY/d(0/H) = 45 ± 19. 

The Pop III star comoving density as a function of redshift, Pi,{z), is given by 



where n{Mh, z) is the comoving number density of dark matter halos of mass Mh at redshift 
z given by the Press & Schechter formalism (Press & Schechter 1974); the integral gives 
the dark matter mass per unit volume contained in dark matter halos with mass greater 
than Mmin{z), where Mmin (computed by Fuller & Couchman (2000)) is a cutoff mass below 
which halos cannot form stars due to the lack of cooling. Qm and fib are the total matter and 
baryon density* in units of the critical density pc — SHq/SttG. The ratio Qb/QM converts the 
integral into baryonic mass which is then turned into stars with an efficiency The latter 
is constrained by the contribution of the redshifted light of Pop III stars to the NIRB. The 
observational data in the NIR (sec e.g. Hauscr & Dwck (2001) and references therein) are 
well fitted by a model in which the contribution from Pop 111 stars is added to the "normal" 
galaxy background light (Salvaterra & Ferrara 2002), where the Pop III contribution is given 



with lu{z, 0) being the specific luminosity of the population at redshift z, computed using 
the spectra for metal- free stars obtained by Schaerer (2002), is the initial mass function 
(IMF), u — uo{l -I- 2;)/(l -I- zq), Te//(^'o, Zo, z) is the effective optical depth at vq of the IGM 
between redshift Zq and z, and dl/dz is the proper line element. The Pop III star formation 
efficiency required to account for the observed NIRB depends essentially only on IMF for 

* We adopt the 'concordance' model values for the cosmological parameters: h = 0.7, CIm = 0.3, Qa = 0.7, = 0.038, 
rjs = 0.9, and F = 0.21, where h is the dimensionless Hubble constant, Hq = 100ft. km Mpc~^ 
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those objects, the redshift z^nd at which the formation of Pop III stars ends being constrained 
to be Zend — 8.8 by the J band data (Salvaterra & Ferrara 2002). In addition to the NIRB, 
an early star formation may contribute to the metal enrichment of the intergalactic medium 
(IGM). Although the enrichment history of the IGM at high redshift is very uncertain, mainly 
depending on the metal escape fraction from galaxies and their actual spatial distribution, it 
sounds reasonable to adopt the mean Lya forest metallicity as a 'reliable' first-order estimate 
of the mean IGM metallicity. Lya forest QSO absorption line experiments indicate a mean 
IGM metalhcity oi\og{ZiGM / Zq) = -2.5±0.5 at ;2 - 3 (Dave et al. 1998, Tclfcr et al. 2002). 
Salvaterra & Ferrara (2002) have shown that this constraint is not exceeded by first stars if 
either [i] their primordial Initial Mass Function is a Salpeter one (0(M) oc M~^-^^), in which 
pair instability supernovae (SN^-y, i.e. progenitor stars with mass in the range 130-260 Mq 
(Heger & Woosley 2002)) are the dominant sources of heavy elements, or [ii] only stars with 
masses in excess of 260 M© are formed, that lock their nucleosynthetic products into a very 
massive black hole (VMBH) remnant (Schneider et al. 2002). In the first case, the NIRB is 
well fitted with — 0.534 and Zend — 8.788. For the second case, although the model seems 
to overproduce the VMBH density to respect of the estimate density of super massive black 
hole (SMBH) (Merritt & Ferrarese 2001), at the moment there is no strong reason to reject 
this model, since the connection between VMBHs and SMBHs is not clear. We adopted here 
0(M) = (5(1000 Mq) and the NIR data are reproduced with = 0.043 and Ze„d = 8.830; 
we will use these Zend values in the following. Note that, as stars more massive than about 
3OOM0 have virtually the same spectrum, the 5-function distribution does not represent an 
unreahstic assumption (Bromm, Kudritzki & Loeb 2001). These two best fit models to the 
NIRB data are shown in Figure 1. Obviously, the NIRB fraction accounted for by PopIII 
stars is dircetly proportional to as a consequence the same linear dependence is found 
for AYp, as it is clear from eq. 3 below and eq. 1. 

The constraints on the IMF obtained from this fit can then be used to estimate the contri- 
bution of Pop III stars to ^He primordial abundance. This is given by 

AYp{4>, z) = fHe{<P) n.(0, Z)/Qb (3) 

where ^^(z) = P*/ Pc] fne is the ^He fraction of the stellar mass ejected via supernova 
explosions and/or through mass- loss in a radiatively-driven wind (although Kudritzki (2002), 
has hinted that winds might not be effective at Z — 0, recent results show that winds can 
take place even in metal-free stars; Lentz, Hauschildt, Aufdenberg & Baron, in preparation); 
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its value depends on the adopted IMF, 0. The ^He produced by SN^^ explosions of Pop III 
stars has been computed by Heger & Woosley (2002). Stars with mass above 260 Mq are 
likely to end up in VMBHs without ^He ejection. On the other hand, according to very 
recent calculations (Marigo, Chiosi & Kudritzki 2002), wind mass-loss is significant only for 
very massive objects (750-1000 Mq). For a 1000 M© non rotating star, 17.33% (0.03%) of 
the initial mass is ejected via mass-loss as helium (oxygen) during the stellar lifetime. A 
small amount of carbon and nitrogen is ejected as well. If rotation is taken into account the 
percentages are 16.97% (He) and 0.15% (O). Pulsational instabilities may not have sufficient 
time to drive appreciable mass-loss in metal free stars (Baraffe, Heger & Woosley 2001). 
We calculate the contribution of Pop HI stars to the ^He abundance according to eq. (3) 
down to redshift Zf,nd at which the formation of Pop HI stars ends. For the Salpeter IMF 
(case [i] above) AYp is quite small, 2 x 10~^ if mass-loss is not considered. This value 
doubles if one accounts for He ejected from 750-1000 Mq stars via wind losses; inclusion 
of rotation does not change these conclusions. In case [ii] there is no contribution from 
SN^^, since 1000 Mq stars fall above the SN^^ progenitor range; as already mentioned, 
mass-loss in a wind is instead quite efficient at ejecting ^He. The resulting contribution is 
AYp = 3.191 X 10~^ (3.258x10"^) in the rotating (non-rotating) case. Moreover, the gas is 
also enriched in oxygen^ to [0/H]= —2.33 (—3.02), leading to a mean IGM metallicity of 
log{Z/ZQ) = —2.52 (—3.15). In Figure 2 we show the ^He data points from the Izotov & 
Thuan sample (Izotov & Thuan 1998) from which we have subtracted the best fit primordial 
value, YpjT- Our results for the non-rotating and rotating cases are also indicated in the 
Figure by the square and triangle points. Both points fall beyond the la error fine (and the 
non- rotating one beyond 2a) . Hence the conclusion is that a measure of the ^He primordial 
abundance extrapolated to zero metallicity from the values of currently sampled [0/H]> 
— 1.6 would include the contribution from the very first generation of stars and therefore 
overestimate the amount of 'truly' primordial, BEN ^He. As a result, a correct determination 
of BBN '^He can only be made by sampling a strictly metal-free gas: any metal-enriched gas 
will already bear the signature of the ^He production of Pop HI stars. A similar conclusion, 
although based on totally different arguments, has been suggested in the past (Mathews, 
Boyd & Fuller 1993). We note that the exact horizontal position of the theoretical points 

t Following common notation, we use square brackets to express values in units of solar abundance, [0/H]=log(0/H) — 
Iog(O/H)0, and use (O/H)© = 5 x 10~*, where (O/H) is the abundance ratio of the two species by number. 
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is affected by the systematic uncertainty introduced by the imphcit assumption made of 
a homogeneous distribution of the stellar nucleosynthetic products; such error is difficult 
to estimate. In general, as observations are pushed to lower metallicity gas, we expect the 
scattering in the data to increase as a result of the inhomogeneous distribution of oxygen, 
which is presumably enhanced close to the most actively star forming regions of the universe. 
We reiterate that the value for AY^ here inferred has to be considered as an upper hmit, 
since a broader IMF in the range 300-1000 Mq leads to a lower contribution. 

The inferred value of AYp gives us a measure of the accuracy with which Yp can be 
measured though current observations of this type. In the most extreme case we can expect 
Yp = YpjT - AYp = 0.2443 - 3.258 x 10"^ = 0.241. Is this low value at odd with BBN 
theory ? It is well known that there is a tension between theoretical BBN predictions and 
observed abundances, in the sense that the latter typically fall short, particularly for ^He and 
^Li. The further decrease of ^He abundance suggested here would aggravate the discrepancy, 
a fact which should lead us to consider even more seriously non-standard BBN models. 
Similar arguments apply if one considers D and "^He data simultaneously: observations of 
high- 2; QSO absorbers seem to suggest a low primordial deuterium abundance, implying high 
nucleon-to-baryon ratio and Yp values. The contribution from Pop III stars will instead tend 
to decrease Yp, a problem which would be exacerbated if we are to believe lower (with respect 
to the IT one) determinations of such quantity, Yp = 0.234±0.003, obtained by other groups 
(Olive, Skillman & Steigman 1997). A lower Yp value would require a slower expansion of the 
radiation-dominated universe. This behavior is predicted by a class of so-called "extended 
quintessence" models, as a result of a decrease of the effective gravitational constant G. 
Obviously, there is a lot that we can learn from BBN element abundance measurements. 
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cation; R. Schneider, E. Skillman and M. Tosi for useful comments and discussions. This 
work was partially supported (AF) by the Research and Training Network 'The Physics of 
the Intergalactic Medium' set up by the European Community under the contract HPRN- 
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Figure 1. 

Figure 1 Caption 

Best fit results of the NIRB data (Salvaterra & Ferrara 2002). Dotted line: Salpeter IMF 
Zend = 8.788, = 0.534. Solid line: = 5(1000 Mq) z^nd = 8.83 = 0.043. The small filled 
circles are the NIRS data (Matsumoto et al. 2000). The open symbols are the DIRBE results: 
squares for Wright (2001), diamonds for Cambresy et al. (2001), and circles for Gorjian et 
al.(2000). The big filled circle is the Kiso star count measurement. The errors are at la and 
for all the data the Kelsall's model (Kelsall et al. 1998) for the zodiacal light is applied. The 
"normal" galaxy contribution is estimated from the Subaru Deep Field (Totani et al. 2001). 
For A > 2.2 jivci a zero contribution from normal galaxies to the NIRB is assumed. 
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Figure 2. 

Figure 2 Caption 

Effect of Pop III stars on the ^He primordial abundance determination. The plot shows the 
^He abundance, Y, to which the estimate by Izotov & Thuan, YpjT (Izotov & Thuan 1998) 
has been subtracted, as a function of the oxygen content of the galaxy. The triangle (square) 
represents the ^He abundance obtained taking into account the helium ejected via mass-loss 
from rotating (non-rotating) Pop III stars. The black dots arc the data of the sample of 
Izotov & Thuan (1998). The solid line is the best fit relation Y - 0/H and the dashed lines 
represent the corresponding la errors (Izotov & Thuan 1998). 
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